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Abstract. The statistical properties of active and passive scalar fields transported 
by the same turbulent flow are investigated. Four examples of active scalar have been 
considered: temperature in thermal convection, magnetic potential in two-dimensional 
magnetohydrodynamics, vorticity in two-dimensional Ekman turbulence and potential 
temperature in surface flows. In the cases of temperature and vorticity, it is found that 
the active scalar behavior is akin to that of its co-evolving passive counterpart. The two 
other cases indicate that this similarity is in fact not generic and differences between 
passive and active fields can be striking: in two-dimensional magnetohydrodynamics 
the magnetic potential performs an inverse cascade while the passive scalar cascades 
toward the small-scales; in surface flows, albeit both perform a direct cascade, the 
potential temperature and the passive scalar have different scaling laws already at 
the level of low-order statistical objects. These dramatic differences are rooted in the 
correlations between the active scalar input and the particle trajectories. The role of 
such correlations in the issue of universality in active scalar transport and the behavior 
of dissipative anomalies is addressed. 
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1. Introduction 

Scalar fields transported by turbulent flows are encountered in many natural phenomena 
and engineering problems, ranging from atmospheric physics to combustion 2 , to 
the transport and amplification of magnetic fields in astrophysical fluids pj. Here, 
we consider the case of advected scalar fields such as temperature, pollutant density, 
chemical or biological species concentration. 

In many cases there is a two-way coupling between the scalar and the flow: the 
transported field can influence the velocity field - this is dubbed active transport. This 
is the case, for example, of the temperature field that acts on velocity through buoyancy 
forces. Conversely, situations where the feedback of the scalar field is negligible and the 
velocity determines the properties of the scalar, but not vice versa are known as passive. 
This ideal case is well approximated by fluorescent dye used in laboratory experiments 
to mark fluid parcels. 

Although active and passive scalars are governed by the same advection-diffusion 
equation, their nature is radically different. Passive scalars belong to the realm of linear 
problems, even though highly nontrivial. Indeed, as a consequence of the statistical 
independence of the forcing and the advecting velocity, the transported field depends 
linearly on the forcing. This property allows a full-fledged theoretical treatment of the 
problem, and has the major consequence that the passive scalar scaling laws are universal 
with respect to the injection mechanism. On the contrary, for active fields the presence 
of the feedback couples the velocity to the transported scalar and makes the problem 
fully nonlinear. In this case, the theoretical tools developed for the study of the passive 
problem may fall short of explaining the behavior of active scalars, and, as of today, the 
understanding of active turbulent transport lags far behind the knowledge accumulated 
for the passive counterpart. This state of the art motivated us to pursue a "case study" 
of turbulent transport of active and passive scalars, using the scaling properties of fields 
evolving in the same turbulent flow as the basic diagnostics for comparison. 

We consider four different systems belonging to the two following general classes of 
problems: (i) active scalars that influence the flow through local forces; (ii) active fields 
functionally related to the velocity. 

The evolution of a scalar belonging to the first class is described in terms of the 
following set of equations: 



where a and c are the active and passive scalar fields, respectively and obey the 
advection-diffusion equations (HJ) and (0). The scalar inputs f a and f c have a 
characteristic lengthscale if, and represent two different realizations of the same 
stochastic process. Were they coincident, no difference between active and passive fields 
would persist. For the sake of simplicity, we use the same molecular diffusivity, k for 



d t a + v ■ Va 
d t c + v ■ Vc 
d t v + v • V« 



k A a + f a , 
k A c + f c , 

- Vp + i/At; + F[a, Va, . . .] , 



(1) 
(2) 
(3) 
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both scalars. The active character of a is embodied by the term F[a, Va, . . .], which 
acts as a forcing for the velocity field in the Navier-Stokes equations. The specific form 
of the F term depends on the physical system under investigation. In the following 
we shall consider two examples from this class: thermal convection j3J where a 
is the temperature field and F = —j3ag is the buoyancy force, and two-dimensional 
magnetohydrodynamics e 6j, where a is the magnetic potential and F = — AaVa is the 
Lorentz force. 

The second class of active scalars is relevant to geophysical flows IHj. In this 
case the dynamics is described in terms of a scalar field obeying the advection-diffusion 
equation supplemented by a functional relation that gives v in terms of a : 



Here, the vector valued kernel T is divergence-free. A well-known instance from such 
class is the two-dimensional Navier-Stokes equation, where the active scalar is the 
vorticity, V x v. Another problem which we shall discuss is the turbulent flow on 
the flat surface of an infinitely high fluid, described by the surface-quasi-geostrophic 
equation [TUj. Here, the active scalar is the fluid density, which is related (e.g. for 
ideal gases) to the so-called potential temperature. 

Some of the results presented in this paper were previously discussed in Refs. QT] 
H21 H3] . Related investigations on active and passive transport may be found in [TJl [TH] 
and [IEIIII!. 

The paper is organized as follows. In Sect. |2J we briefly review some results 
about passive scalar transport in turbulent flows. Particular emphasis is put on the 
Lagrangian description of scalar transport, pointing out the results which hold for 
active scalars as well. In Sect. 01 the statistics of the (active) temperature and a 
passive scalar field in 2c? convection is discussed. Sect. HI is devoted to the analysis 
of two-dimensional magnetohydrodynamics, with a thorough discussion of the problem 
of dissipative anomaly in scalar transport. In Sect.Elthe evolution of passive and active 
fields in Ekman-Navier-Stokes turbulence is discussed. In Sect.|H]we study the properties 
of turbulence on fluid surfaces under the quasi-geostrophic approximation. In the last 
section we summarize the main findings and delineate the perspectives for active scalar 
transport. 

2. Passive scalar turbulence 

2.1. Eulerian approach 

The dynamics of passive scalars is governed by the advection diffusion equation @. 
In order to describe the general properties of passive scalar evolution in turbulent 
incompressible velocity fields, we assume that the velocity field v is scale invariant 
and rough. In other words, the spatial increments 8 r v = (v(x + r, t) — v(x, t)) ■ r/r 
depend on the separation fractional power, i.e. 5 r v ~ r with h < 1 (e.g. 




(4) 
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h = 1/3 in Kolmogorov's 1941 turbulence fTS]). Being interested in the statistically 
steady properties of the field, we introduce a source of scalar fluctuations f c . In the 
following we take for convenience a random, Gaussian, statistically homogeneous and 
isotropic forcing with zero mean and correlation function 



The correlation function of the forcing Tirjif) is roughly constant at scales smaller 
than if, which is assumed to be within the scaling range of v, and drops rapidly to zero 



The phenomenology of passive scalar turbulence may be summarized as follows. 
Scalar fluctuations injected at the scale if are transferred toward the small scales with 
a constant flux, down to the dissipative scale id- There, the molecular diffusion absorbs 
the incoming flux and ensures the equilibrium between the input and the dissipation. 
The fluctuations are thus maintained in a statistically steady state which is characterized 
by two major properties. First, the scalar dissipation is asymptotically independent of 
the molecular diffusivity k, attaining a finite nonzero limiting value for k — > 0. This 
singular behavior of the dissipation is also known as dissipative anomaly. Second, in the 
scaling range Id <C r <C if the scalar statistics is intermittent. This amounts to saying 
that the small-scale statistics is characterized by the alternance of strong, rare events 
where the scalar increments, 5 r c = c(x + r,t) — c(x,t) are much larger than their typical 
value, c rms , and long quiescent phases where 5 r c <C c rms . Intermittency is reflected by 
the scaling behavior of the structure functions, i.e. the moments of the scalar increments 



The scaling exponents Cat are s& id to be anomalous, when they deviate from the 
dimensional expectation = N(l — h)/2. The equality = (f}" 1 holds possibly 
only for N = 2 4j, whereas for N > 2 the deviations become more and more severe. 
The asymptotic behavior at large orders N corresponds to the saturation — > 
[2T] . The saturation is related to the presence of sharp "fronts" in the scalar field. The 
exponents are universal with respect to the details of the energy injection statistics. 
The forcing only affects the numerical prefactors appearing in the structure functions. 

The anomalous scaling <^ m ^ signals the breakdown of scale invar iance, as 
confirmed by the explicit appearance of if in (jHJ), even at scales r -C if- Indeed, 
anomalous scaling of the moments of scalar increments is equivalent to state that 
the probability density functions (pdf's) of 5 r c at different r's cannot be collapsed 
by rescaling them with a unique power law in r. Even though the specific values of 
the exponents depend on the details of the flow (statistics, time correlation, and 
roughness exponent h) , intermittency and the breaking of scale invariance in the scalar 
statistics are generic features of passive scalar turbulence. 

The physical mechanism leading to anomalous scaling has been recently understood 
in the framework of the Kraichnan model of scalar advection [22123 (see Refs. P4*1I25| I2B] 
and Ref. [2Zj for an exhaustive review on the subject). In this model the advecting 



{Mxi, t)f c (x 2 , t')> = S(t - t>)F(\ Xl - X 2 \/if) . 



(5) 



for r > if. 




dim 



(6) 
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flow v is random, Gaussian, self-similar and 5-correlated in time. Under these special 
conditions, there exists a closed set of linear equations for multi-point correlation 
functions. The anomalous exponents are the scaling exponents of the homogeneous 
solutions (the so-called zero modes) of those equations. Since homogeneous solutions 
do not depend by definition on the scalar input, their scaling exponents are universal, 
and cannot be inferred from dimensional arguments. The concept of zero mode can be 
extended to passive scalar turbulence in generic velocity fields 

The properties of passive scalars described up to now are in the language of fields 
-the Eulerian description. It is now interesting to adopt a different but equivalent 
viewpoint in terms of particle trajectories -the Lagrangian description. 

2.2. Lagrangian description 

The basic idea of the Lagrangian approach is to solve Eq. (J2J) by the method of 
characteristics. Let us denote p(s; x, t) (henceforth whenever there is no ambiguity 
we indicate it as p(s)) as the trajectory of a fluid particle landing at point x at time t. 
The path p(s) is the solution of the stochastic differential equation 
dp{s) 



v( P (s),s) + V2kw{s), p{t) = x, (7) 



ds 

where i,j = l,...,d(d being the space dimensionality), and w(s) is a Wiener process 
(the derivative of a Brownian motion), i.e. Wi are Gaussian variables of zero mean and 
correlation {wi(s)ijbj(s')) = Sij5(s — s'). Along the path p(s) Eq. (J2J) reduces to 

JL±± = f e (p(s),s), (8) 

which is easily solved as <p w {t) = J" * ds / c (p(s), s). For the sake of simplicity, we assumed 
0™(O) = 0. We indicated with (fi w the solution obtained along the path p obtained for 
a specific realization of the process w. The passive scalar field, c(x,t), is recovered by 
averaging over all the realizations of w, i.e. along all the Lagrangian paths ending in x 
at time t 



c(x,t) = (<p w (t)) w =( / dsf c (p(s),s)) . (9) 







The statistic of the trajectories is summarized in the particle propagator 
P(y,s\x,t) = (5(y — p(s;x,t)) w , that is the probability of finding a particle at point 
y and time s < t, provided it is in x at time t. According to the theory of stochastic 
processes (HO], P(y,s\x,t) obeys the Kolmogorov equations: 

- d s P(y,s\x,t) - V y ■ [v(y,s)P(y,s\x,t)] = KA y P(y,s\x,t) , (10) 

d t P(y,s\x,t) + 'V x -[v(x,t)P(y,s\x,t)] =K,A x P(y, s\x, t) , (11) 

with initial condition P(y,t\x,t) = 5(x — y). The unusual minus signs in the l.h.s. of 
fT0|) are due to the fact that particles move backward in time. The solution of (J2J) can 
be written in terms of the propagator: 

c{x,t) = ds J dy f c {y,s)P(y,s\x,t) , (12) 
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as it can be directly checked by inserting (f*H2j) in (J2J) and using (fTTj). 

At variance with smooth velocities (i.e. Lipschitz continuous, 5 r v ~ r) where for 
k — > particle trajectories are unique, for velocity fields rough and incompressible 
the particle propagator does not collapse onto a single trajectory in the limit k — > 0. 
Lagrangian paths are not unique and initially coincident particles separate in a finite 
time. This property is at the root of the dissipative anomaly. 

For active and passive scalars evolving in the same flow, the Lagrangian paths p(s) 
are the same, as well as the propagator P(y, s\x,t). The difference between a and c is 
that, since the active scalar enters the dynamics of v (013), the Lagrangian trajectories 
are functionally related to the active scalar forcing f a , but are independent of the passive 
source term f c . This is the Lagrangian counterpart of the linearity of the passive scalar 
problem, which does not hold for the more complicated case of active scalars. 

One of the advantages of the Lagrangian description is that it makes transparent 
the physics of transport processes. For instance, let us consider the 2-points correlation 
function for passive scalars, C c (xi, x 2 \ t) = (c(x\, t)c(x2, t)). This offers the possibility 
of an intuitive interpretation of the energy cascade phenomenology and gives insights 
into measurable statistical objects such as the scalar spectrum E c (k), that is the Fourier 
transform of C 2 . From Q and averaging over f c and v one obtains 

Cl{x x ,x 2 \t) = (/ dsi / ds 2 {f c (p(s 1 ;x 1 ,t))f c (p(s 2 ;x 2 , t))) f ) . (13) 

WO JO I W v 

Introducing the velocity-averaged two-particle propagator (P 2 (yi, y 2 , s\xi, x 2 , t)) v , 
which evolves according to the straightforward generalization of (jl()| lll )) to two particles, 
and using ©, Eq. (fT^jl reduces to 

C 2 (x h x 2 ; t) = J ds J J (P 2 (yi,y2,s\x 1 ,x 2 ,t)) v J r (\y 1 -y 2 \/e f )dy 1 dy 2 , (14) 

which has a clear physical interpretation. Since J-(x) is vanishingly small for i>1, the 
correlation C 2 is essentially the average time spent by a particle pair with a separation 
r = \x\ — x 2 \ below the forcing scale if. Due to the explosive separation of particles, 
this time has a finite limit for r — > 0, which yields the leading contribution to C 2 . The 
subleading behavior is uncovered by the second-order structure function 

S c 2 (r) = ((c( Xl ,t) - c(x 2 , t)) 2 ) = 2(C 2 C (0) - C 2 c (r)) , (15) 

which is roughly the average time T~i f (r) taken by two coinciding particles to reach a 
separation r. For a Kolmogorov 1941 turbulent flow (h = 1/3) one has S 2 (r) ~ r 2 / 3 , i.e. 
E c {k) ~ A;~ 5 / 3 - the Oboukov-Corrsin dimensional expectation jlj. 

The Lagrangian description can be extended also to higher-order statistics as 
multipoint correlation functions C^(xi, . . . , x N ) = (c(xi, t) . . . c(x N , t)). However, 
when many points come into play, their geometrical arrangement becomes crucial. 
Dimensional arguments, which are based on the size of the configuration but forcibly 
neglect the "angular" information, fall short of capturing the observed behavior for 
multi-point observables. A detailed discussion of their properties is beyond the scope 
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of this brief review. In the following we just summarize the main concepts, referring to 
Refs. jHH HI] for further reading. 

Expanding the power in the definition of the structure functions S%(r), it is 
immediate to express them as a linear combinations of iV-point correlation functions 
(see, e.g., (fTHjl for S'f(r)). Therefore, the latter must contain a contribution, denoted 
as Z^(xi, . . . ,x N ), that carries the anomalous scale dependence: Z^(Xxi, . . . , Xx N ) = 
X^Z^i^Xi, . . . , x N ) |2U 123 EH]- From a Lagrangian viewpoint, the function Z^ has a 
special property that distinguishes it from a generic scaling function. The remarkable 
result is that |3H I2E] 

j t (Z%) c = 0, (16) 

where the derivative 4 is performed along the trajectories of N particles advected by 
the flow, and the average is over the ensemble of all trajectories. In other terms, Z C N 
is statistically preserved by the flow [2H1 E2j Universality of scaling exponents is then 
just a byproduct of the definition of statistically preserved structures: since f c does not 
appear in Eq. (jlfij) . the properties of zero modes are insensitive to the choice of the 
forcing. 



2.3. Dissipative anomaly 

In spite of the continuous injection of scalar through the pumping f c , the second-order 
moment (c 2 (x, t)) does not grow indefinitely even in the limit k — > 0. This is due to the 
existence of a finite nonzero limit of the scalar dissipation e c = k\ Vc| 2 - the dissipative 
anomaly. 

In order to understand how (c 2 (x,t)) achieves a finite value independent of the 
diffusivity coefficient, we adopt the Lagrangian viewpoint. From Eq. (jlHj) we have: 

(c 2 (x,t))=Q^ jf ds 1 ds 2 f c (p(s 1 ;x,t))f c (p(s 2 ;x,t))^ = / dsf c {p{s; x, \ ,(17) 

where the brackets indicate the average over the scalar forcing, the velocity field and 
the noise. 

Looking naively at (|17|) one might expect that for a large class of random forcing 
of zero mean the r.h.s of the above expression would grow linearly with t. For instance, 
when the forcing is Gaussian and 5-correlated in time, one could argue that (fTTj) is 
essentially a sum of independent variables and by central limit theorem arguments 
conclude that (c 2 ) oc t. This conclusion would be correct if in the limit k — > all 
trajectories collapse onto a unique Lagrangian path. This turns out to be the case for 
strongly compressible flows but not in general. For compressible flows, energy indeed 
grows linearly in time and the advected scalar performs an inverse cascade process 
On the contrary, in rough incompressible flows coinciding particles typically 
separate in a finite time, giving rise to multiple paths. As a consequence, a self-averaging 
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process takes place in (|17|) and this prevents the indefinite growth of the energy. This 
is evident upon rewriting (j!7|) as 

(c 2 {x, t)) = J ds j j (Pad/j, y 2 , s\x, x, t)) v F{\ yi - y 2 \/£ f ) d yi dy 2 . (18) 

The time integral is cut off at \t — s\ ^> % f , that is for times larger than the (finite) 
time needed by two coinciding particles to separate by a distance larger than the forcing 
correlation length if. This is the mechanism leading to a finite dissipation of energy. To 
summarize, the incompressibility and the roughness of the flow result in the dissipative 
anomaly through the explosive separation of particle paths. Further discussion on the 
role of dissipative anomaly in passive scalar turbulence can be found in [23 EHJ ES] ■ 

3. Two-dimensional turbulent convection 

An interesting problem in the context of turbulent transport is the advection of 
inhomogeneous temperature fields in a gravitational field. Temperature fluctuations 
induce density fluctuations that in turn, via buoyancy forces, affect the velocity field: 
hence the temperature field is an active scalar |UEj- Here, we consider two-dimensional 
convection, which is also of experimental interest in Hele-Shaw flows jSE]. As an 
additional asset, the two-dimensional problem is better suited for the study of scaling 
properties since it allows to reach higher resolution and larger statistics. 
Two-dimensional convection is described by the following equations 

d t a + v-Va = kA a + f a , (19) 
d t v + v ■ Vv = — Vp + vAv — Pag — av , (20) 

where a is the field of temperature fluctuations. The second equation is the two- 
dimensional Navier-Stokes equation where v is forced by the buoyancy term —f3ga 
in the Boussinesq approximation g — gy is the gravitational acceleration and f3 
the thermal expansion coefficient. Kinetic energy is removed at the large scales by the 
friction term, —av. The friction is physically due to the drag experienced by a thin 
(quasi 2d) layer of fluid with the walls or air [7J |37] ; ct is related to the thickness of the 
fluid layer. A passive scalar, c, evolving according to Eq. (J2J) in the same flow has been 
considered as well, for comparison. 

Before entering in the active/passive scalar issue, let us briefly recall the 
phenomenology of 2D turbulent convection (for the 3D case see, e.g., |U El)- The 
balance of buoyancy and inertial terms in (|2*U|) introduces the Bolgiano lengthscale 
C.B- \\- At small scales, r £b, the inertial term is larger than buoyancy forces 
and the temperature is basically a passive scalar. At large scales, r 3> £ B , buoyancy 
dominates and affects the velocity, which in two dimensions performs an inverse 
energy cascade. However, at variance with the usual 2D Navier-Stokes turbulence, 
the kinetic energy input rate e depends here on the scale. Dimensional arguments yield 
s(r) = (3g ■ (v(x + r,t)a(x,t)) ~ r 4//5 , the Bolgiano scaling for the velocity structure 
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Figure 1. Scaling exponents of temperature, (fa, and velocity, (fa. The straight lines 
are the dimensional predictions, N/5 for temperature, 3iV/5 for velocity. Details on 
the numerics may be found in [U]. Notice that at orders larger than N = 8 the 
temperature exponents saturate to a constant value ~ 0.8. 



functions 

S;(r)~(#f/ 3 ~/«, G = 3iV/5, (21) 

and for temperature 

S£(r)~r<», G = iV/5. (22) 

No intermittency corrections are reported for the velocity, whereas the temperature field 
is strongly intermittent (see Fig. [T] and Refs. jHEllII])- 

Summarizing, the temperature fluctuations are injected at scales, ~ if, pump 
kinetic energy through the buoyancy term, and a non-intermittent velocity inverse 
cascade establishes with 5 r v ~ r 3 / 5 . The presence of friction stabilizes the system 
inducing a statistically steady state. The active and the passive scalars are therefore 
transported by a self-similar, incompressible and rough flow. An intermittent cascade of 
fluctuations with anomalous, universal scaling exponents, (fa is observed for the passive 
scalar. 

Our aim is to compare the statistical properties of the temperature and the passive 
scalar fields. To this purpose, in Ref. [TT], Eqs. (fTT?j) and (J2UJ) have been integrated with 
f a and f c , chosen as two independent realizations of a stochastic, isotropic, homogeneous 
and Gaussian process of zero mean and correlation: 

(f t {x,t)f 3 {x',t')) = 5 tJ 5{t-t')F{\x-x'\/if), (23) 

where T{rjlf) = exp(— r 2 /(2£j? )) decreases rapidly as r > if. The labels are i,j = a,c. 
The results of this numerical study clearly confirm that temperature scaling exponents 
are anomalous (Fig. [TJ), and coincide with those of the passive field: (fa = (fa (Fig. |2k>). 

In this system there is saturation of intermittency, i.e. for large iV the scaling 
exponents saturate to a constant (^ c ~ 0.8 (see Fig. QJ). This phenomenon, well known 
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Figure 2. (a) Local scaling exponents of temperature (x) and concentration (A) 
fluctuations, Cjv C ( r ) = rflnS'^/ c (r)/dlnr. Temperature and concentration are driven 
by independent Gaussian random forcing with correlation function as Eq. (J2HJ). (b) 
Temperature structure functions, S^(r), for N — 2, 4, 6, as a function of the separation 
r. The two sets of curves are generated by using the random ( x ) and mean gradient (+) 
forcing. Note the parallelism within the inertial range. Curves have been multiplied 
for a numerical factor for visualization purposes. 

for passive scalar s 21,, is physically related to the presence of abrupt changes in the 
spatial structure in the scalar field ("fronts"). In the temperature field these quasi- 
discontinuities correspond to the boundaries between hot rising and cold descending 
patches of fluid j^HJ. It is worth mentioning that saturation has been experimentally 
observed both for passive scalars [THl EOj and temperature fields [HE] • Clear evidences 
of saturation have recently been obtained also in the convective atmospheric boundary 
layer exploiting the large-eddy simulation technique [4*U] . 

These findings point to the conclusion that the temperature and the passive scalar 
have the same scaling laws. It remains to be ascertained whether the temperature scaling 
exponents are universal with respect to the forcing. To this aim, a set of simulations has 
been performed in Ref. jSHlEj, with a forcing that mimics the effect of a superimposed 
mean gradient on the transported temperature field 

f a {x,t) = yg • v(x,t) . (24) 

Remarkably, the results show that the scaling exponents of the temperature field do not 
depend on the injection mechanism (Fig. Eb) suggesting universality [TT| 138]. Another 
outcome of this investigation is that the velocity field statistics itself is universal with 
respect to the injection mechanism of the temperature field. Indeed, v displays a close- 
to-Gaussian and non-intermittent statistics with both forcing (f2*HJ) and (|2"4"|) [TTj . This 
is most likely a consequence of the observed universal Gaussian behavior of the inverse 
energy cascade in two-dimensional Navier-Stokes turbulence |41[ l4"2~] . Indeed, velocity 
fluctuations in two-dimensional convection also arise from an inverse cascade process 
driven by buoyancy forces. 

So far all the numerical evidences converge to the following global picture of scaling 
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and universality in two-dimensional turbulent convection. Velocity statistics is strongly 
universal with respect to the temperature external driving, f a . Temperature statistics 
shows anomalous scaling exponents that are universal and coincide with those of a 
passive scalar evolving in the same flow. It is worth noticing that similar findings have 
been obtained in the context of simplified shell models for turbulent convection [T4"|ll5j. 

The observed universality of the temperature scaling exponents suggests that a 
mechanism similar to that of passive scalars may be at work, i.e. that statistically 
preserved structures might exist also for the (active) temperature. Pursuing this line of 
thought, one may be tempted to define them through the property c = as for 

passive scalars (see Eq. (fTHjl ). However, statistically preserved structures are determined 
by the statistics of particle trajectories which, through the feedback of a on v, depend 
on f a . Therefore, the above definition does not automatically imply the universality of 
Zfj, because Lagrangian paths depend on f a . Nonetheless, the observed universality of 
the statistics of v is sufficient to guarantee the universality of the trajectories statistics, 
leading to the conclusion that if Zfj exists it might be universal. Since Z^ are also 
defined by the Lagrangian statistics that is the same for a and c, one might further 
conjecture that Z^ = Z%. This would explain the equality of scaling exponents Cn = Cn- 

It has to be remarked that this picture is not likely to be generic. Two crucial points 
are needed to have the equality between active and passive scalar exponents: (i) the 
velocity statistics should be universal; (ii) the correlation between f a and the particle 
paths should be negligible. As we shall see in the following, those two requirements are 
not generally met. 

4. Two-dimensional magnetohydrodynamics 

4-1. Direct and inverse cascades 

Magnetohydrodynamics (MHD) models are extensively used in the study of magnetic 
fusion devices, industrial processing plasmas, and ionospheric/astrophysical plasmas 
MHD is the extension of hydrodynamics to conductive fluids, including the effects of 
electromagnetic fields. When the magnetic field, 6, has a strong large-scale component 
in one direction, the dynamics is adequately described by the two-dimensional MHD 
equations jHj- Since the magnetic field b(x,t) is solenoidal, in 2d it can be represented 
in terms of the magnetic scalar potential, a(x,t), i.e. b = — V^a = (— 82a, dia). The 
magnetic potential evolves according to the advection-diffusion equation 

d t a + v ■ Va = kA a + f a , (25) 

and will be our active scalar throughout this section. The advecting velocity field is 
driven by the Lorentz force, (V x b) x b = — AaVa, so that the Navier-Stokes equation 
becomes 

d t v + v ■ Vv = - Vp + uAv - AaVa . (26) 

The question is whether the picture drawn for the temperature field in 2c? convection 
applies to the magnetic potential as well. 
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Figure 3. Power spectra of active (red) and passive (blue) scalar variances E a (k) = 
irk\a(k, t)\ 2 and E c (k) = 7rfc|c(fc, t)\ 2 . In the lower left corner, the fluxes of scalar 
variance II a , c out of wavenumber k. Negative values indicate an inverse cascade. In 
the upper right corner, the total scalar variance e a c (t) = J E a . c (k,t)dk. The active 
variance e a (t) grows linearly in time whereas e c (t) fluctuates around a finite value (see 
text). The rate of active to passive scalar dissipation is e a /e c ~ 0.005. All fields are 
set to zero at t = 0, and time is defined in units of eddy-turnover time T = lf/v Tms 
where If — 2-rt/kf. At k < kf we observe power-law behaviors E a (k) ~ fc- 2 -°± - 1 an d 
E c {k) ~ k°- 7±0A , while at k > k f we find E a {k) ~ fc-3.6±o.i and E ^ _ fc-i.4±o.i_ 




Eqs. (j25|) and (|2fi|) have two quadratic invariants in the inviscid and unforced limit, 
namely the total energy | J (v 2 +b 2 )dx and the mean square magnetic potential \ J a 2 dx. 
Using standard quasi-equilibrium arguments an inverse cascade of magnetic potential 
is expected to take place in the forced and dissipated case jlSJ- This expectation has been 
confirmed in numerical experiments j^j. Let us now compare the magnetic potential 
with a passive scalar evolving in the same flow. 

We performed a high-resolution (4096 2 collocation points) direct numerical 
simulations of Eqs. (|25j) - (j26|) along with a passive scalar (J2J). The scalar forcing terms f a 
and f c are homogeneous independent Gaussian processes with zero mean and correlation 

(/*(*, O/iC^O) = j2^) 5 v 5 ( k + k '^ k - k ^ 5{t - o ( 2? ) 

where i,j = a,c. The injection length scale If ~ 27r//c/ has been chosen roughly in the 
middle of the available range of scales. F is the rate of scalar variance input. 

In Fig. El we summarize the spectral properties of the two scalars. The emerging 
picture is as follows. While a undergoes an inverse cascade process, c cascades downscale. 
This striking difference is reflected in the behavior of the dissipation. The active 
scalar dissipation, e a = ft|Va| 2 , vanishes in the limit k —>■ - no dissipative anomaly 
for the field a. Consequently, the squared magnetic potential grows linearly in time 
£a(t) — | J a 2 dx fa ^F t. On the contrary, for the passive scalar, a dissipative anomaly 
is present and e c = k|Vc| 2 equals the input holding c in a statistically stationary 
state. (See inset in Fig. EJ). 
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Figure 4. Velocity power spectrum E v (k) — 7rfc|u(fc, t)\ 2 . For k < kf, in agreement 
with previous simulations , we observe E v (k) ~ fc 1 / 3 (dashed lines) which deviates 
from the dimensional prediction fc -1 / 3 . In the range k > kf a scaling close to fc~ 5 / 3 
is observed, indicating that v is rough both in the inverse and in the direct cascade 
range. 



The velocity field is rough (as confirmed by its spectrum, see Fig. 0J) and 
incompressible, therefore particle paths are not unique and explosively separate. This 
entails the dissipative anomaly for passive scalars. On the contrary, the dissipative 
anomaly is absent for the magnetic potential, in spite of the fact that the trajectories 
are the same - the advecting velocity is the same. How can these two seemingly 
contradictory statements be reconciled ? 

4-2. Dissipative anomaly and particle paths 

The solution of the riddle resides in the relationship between the Lagrangian trajectories 
and the active scalar input f a . These two quantities are bridged by the Lorentz force 
appearing in (j2fij) . To study the correlations between forcing and particle paths we need 
to compute the evolution of the particle propagator. The relevant observables are the 
time sequences of 4> a ,c(s) = J dy f a , c (y, s )P(y, s\x,t). Indeed the equivalent of Eq. (fT2|) 
can be written for the active scalar as 



The main difficulty encountered here is that P evolves backward in time according 
to (|10j) and the condition P(y,t\x,t) = 5(y — x) is set at the final time t. Conversely, 
the initial conditions on velocity and scalar fields are set at the initial time. The solution 
of this mixed initial/final value problem is a non-trivial numerical task. To this aim, we 
devised a fast and low-memory demanding algorithm to integrate Eqs. ()25|) - (|26J) and 
fllOjl with the appropriate initial/final conditions. The details are given in Ref. |44j . 

The typical evolution of the propagator is shown in the central column of Fig. |S| 
From its evolution we reconstructed the time sequences of the forcing contributions 
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Figure 5. Time runs from bottom to top. First column: time evolution of the active 
scalar field resulting from the numerical integration of Eqs. 12">l) and l|2t)|) . Second 
column: backward evolution of the particle propagator according to Eq. (|10|l . Third 
column: time evolution of the passive scalar field in the same flow. 

4>a,c(s) which, integrated over s, give the amplitude of the scalar fields according to 
([E} and (J2HJ) . The time series of (j) a (s) and <j) c {s) are markedly different (Fig. EJ, 
the former being strongly skewed toward positive values at all times. This signals that 
trajectories preferentially select regions where f a has a positive sign, summing up forcing 
contributions to generate a typical variance of a of the order F t. On the contrary, the 
passive scalar sequence displays the usual features: f c is independent of P and their 
product can be positive or negative with equal probability on distant trajectories. This 
ensures that the time integral in Eq. (fH2j) averages out to zero for \s — t\ > % f (see 
Sect. E2) and yields c 2 ~ F T tf . 

As shown in the lower panel of Fig. EJ the effect of correlations between forcing 
and propagator is even more striking comparing J"^ ds'</> a (s') with ds'0 c (s'). The 
conspicuous difference has to be related to a strong spatial correlation between P and 
a, as can be inferred from 

J <j) a (s')ds' = j dya(y,s)P(y,s\x : t), (29) 

which can be derived from ()25)1 and An equivalent relation holds for c as well. 

The comparison of the first and the second column of Fig. 03 highlights the role of 
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Figure 6. Top: <j) a . c (s) = f dy f a . c (y, s)P(y, s\x,t). The two graphs have the same 
scale on the vertical axis. Here, t — 32. Bottom: time integrals Jq (f> a (s') ds' (upper 
curve) and f^<p c {s')ds' (lower curve). Note the different scale on the vertical axis. 
Recall that L <j> a (s') ds' — a(x, t) (and similarly for c). 



spatial correlations: the distribution of particles follows the distribution of the active 
scalar. This amounts to saying that large-scale scalar structures are built out of smaller 
ones that coalesce together j3H]- This has to be contrasted with the absence of large- 
scale correlations between the propagator and the passive scalar field (second and third 
column of Fig. [H]) . 

Let us now clarify the mechanism for the absence of dissipative anomaly. Consider 
the squared active field a 2 . It can be expressed in two equivalent ways. On one hand, 
it can be written as the square of On the other hand, multiplying (J2*5j) by 2a one 

obtains the equation 

d t a 2 + v ■ Va 2 = nAa 2 + 2a f a - 2e a . (30) 

Exploiting the absence of dissipative anomaly, e a = 0, Eq. (|5U|) reduces to a transport 
equation that can be solved in terms of particle trajectories. The comparison of the two 
expressions yields jlO] 

J dsj^ ds'jj f a (y, s)f a (y', s')P(y, s\x, t)P(y', s'\x, t) = 

J dsj^ ds'JJ f a (y, s)f a (y', s')P(y', s'; y, s\x, t) , (31) 

where P(y',s';y,s\x,t) = P(y, s\x,t)P(y', s'\y, s) denotes the probability that a 
trajectory ending in (x,t) were in (y,s) and (y',s f ). Integration over y and y' is 
implied. Eq. (|3*T|) amounts to saying that 



f a (p(s),s)ds 



2 



fa(p(s),s)ds 



2\ 



(32) 
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meaning that J * f a (p(s), s)ds is a non-random variable over the ensemble of trajectories. 
The above procedure can be generalized to show that (f f a (p(s),s)ds)^ } = 
([f*f a (p(s), S )d S ] N } w . 

In plain words, the absence of the dissipative anomaly is equivalent to the 
property that along any of the infinite trajectories p(s) ending in (x, t) the quantity 
f* f a (p(s), s)ds is exactly the same, and equals a(x,t). Therefore, a single trajectory 
suffices to obtain the value of a(x,t), contrary to the passive case where different 
trajectories contribute disparate values of J* f c (p(s), s)ds, with a typical spread e c t, 
and only the average over all trajectories yields the correct value of c(x,t). In the 
unforced case particles move along isoscalar lines: this is how non-uniqueness and 
explosive separation of trajectories is reconciled with the absence of dissipative anomaly. 

Inverse cascades appear also for passive scalars in compressible flows jHI]- There, the 
ensemble of the trajectories collapses onto a unique path, fulfilling in the simplest way 
the constraint ()31|1. In MHD the constraint is satisfied thanks to the subtle correlation 
between forcing and trajectories peculiar to the active case. 

Magnetohydrodynamics in two dimension represents an "extreme" example of the 
effect of correlations among Lagrangian paths and the active scalar input. The property 
that all trajectories ending in the same point should contribute the same value of the 
input poses a global constraint over the possible paths. 

Before discussing the statistical properties of a, on the basis of the previous 
discussion it is instructive to reconsider the concept of dissipative anomaly in general 
scalar turbulence. 



4-3. Dissipative anomaly revisited 

In this subsection we give an alternative interpretation of dissipative anomaly. To this 
aim, let us denote with Q(x,t) a generic scalar field, regardless of its passive or active 
character. The scalar evolves according to the transport equation 

d t Q + v • V9 = kAQ + f . (33) 

We can formally solve Eq. (|3*3*|) by the method of characteristics, i.e. in terms of the 
stochastic ordinary differential equations 

dp(t) 

dt 
d#(t) 



v(p(t),t) + V2KW, (34) 
fe(p(t),t). (35) 



dt 

At variance with the procedure adopted in Sect. 12.21 here we are not conditioning a 
priori the paths to their final positions. The Eulerian value of the field is recovered once 
the average over the all paths, p, landing in (x,t) is performed, i.e. Q(x,t) = {&(t)) w . 
Recall that if the flow is non-Lipschitz continuous, such paths do not collapse onto a 
single one also for k — > 0. 

We can now define P(x, t\x , #o, 0) as the probability that a path which started 
in x at time with $ = Q(xq, 0) and arrives in x at time t carrying a scalar value 
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Figure 7. J" dec cr 2 (a;, (boxes) and J dx a^(x,t) (circles) vs time. We have integrated 
Eqs. I|25[) and 10 for two different values of the diffusivity k — 0.003 (filled symbols) 
at resolution 512 2 and n — 0.001 (open symbols) at resolution 1024 2 . The straight 
lines indicates the growth laws 2e a t and 2e c t for the two values of k. Note that e c 
does not depend on n consistently with the presence of dissipative anomaly, while 
e a decreases as k decreases. The variance has been evaluated averaging over 10 6 
Lagrangian paths evolving according to (|34fl . The Lagrangian scalar values t? and 
i9 c (where a(x,t) — f $ a 'P(x,i!) a ,t)d'd a , the equivalent relation holds between c and 
$ c ) have been computed integrating (|3*5)) both for the active f a and passive f c forcings 
along each paths. The forcings are chosen as in (|27[1 The initial Eulcrian a(x, 0), c(x, 0) 
and Lagrangian d a (0), $ c (0) fields have been set to zero. Times is measured in eddy 
turnover times. 



Notice that the conditioning is now on the initial value, and P evolves according to the 
Kolmogorov equation: 

d t P + v ■ V X P + /eVtfP = kAP (36) 

with initial condition P(x, t\x , $ , 0) = S(x — x )5($ — i? ) an d $o — ®(x, 0). 

Integrating over the initial conditions we define now the probability density 
V(x,'d,t) = J P(x, t\x , Q(x , 0), 0)daJo- V still obeys to (J3T?)) with initial condition 
V(x, 0) =5('d — Q(x , 0)), and represents the probability that a path arrives in (x, t) 
carrying a scalar value i?(t) = L* fo(p(s), s)ds. 

Let us now look at the variance of the distribution of such values, i.e. cx e (:E,t) = 
J d 2 Vd{) — (J dVdd) 2 . From Eq. (f3T)J) it is easy to derive the following equation 

d t al(x,t)+v ■V x al(x,t) = KAal(x,t)+2e @ (x,t) (37) 

where e («,t) = k\'V x J -^^(a?, i9, t)d^| 2 . In the Eulerian frame e e (x, t) = k\VQ{x, t)\ 2 , 
i.e. the local dissipation field. 

Integrating over x, we end up with 

j t J a 2 e (x, t)dx = 2 J e e (x, t)dx = 2e , (38) 

where e = (ac| VQ(x, i) \ 2 ) is the average dissipation rate of (0 2 )/2. Therefore, if 6 2 
cascades toward the small scales with a finite (even in the limit k, —* 0) dissipation, ee, 
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Figure 8. (a) The surface is the instantaneous Eulerian magnetic potential field 
a(x, t). The dots represent the positions of particles in the (x, $ a )-space at time t. (b) 
The same for the passive scalar, c. The time t corresponds to the largest time in Fig. EI 
For visualization purposes only the portion [0,7r/4] x [0, 7r/4] is displayed. In (a) the 
cloud of dots closely follows the the magnetic potential surface, while in (b) they are 
considerably more dispersed (notice the differences in the values of the $-axis). 



the variance of the distribution of the values of = J ds/e(p(s), s) will grow linearly 
in time. Conversely, for an inverse cascade of 2 in the absence of dissipative anomaly 
/ dx<jQ(x,t) = 0, corresponding to a singular distribution V(x,$,t) = 5({} — Q(x,t)). 

In Fig. [7|we show the time evolution of J dx a 2 c (x, t) in the case of 2c? MHD which 
confirms the above findings. 

The fact that, for inverse cascading scalars, the probability density V collapses 
onto a 5-function in the limit of vanishing diffusivity amounts to saying that particles 
in the (a;, $) space do not disperse in the ^-direction but move remaining attached to 
the surface i? = Q(x,t). (see Fig. EK) - This is related to the strong spatial correlations 
between a and the Lagrangian propagator observed in Fig. El On the contrary, for a 
direct cascade of scalar such correlations do not exist and dissipation takes place because 
of dispersion in the ^-direction (see Fig. |Ht>)- 



4-4- Eulerian statistics 

4-4-1 ■ Single-point statistics The statistical properties of the magnetic potential a 
and the passive scalar c are markedly different already at the level of the single- 
point statistics. The pdf of a is Gaussian, with zero mean and variance F$t (Fig. EJ). 
Conversely, the pdf of c is stationary and supergaussian (see Fig. EJ), as it generically 
happens for passive fields sustained by a Gaussian forcing in rough flows |27j . 

The Gaussianity of the pdf of a is a straightforward consequence of the vanishing 
of active scalar dissipation. This is simply derived by multiplying Eq. (}25|) for 2na 2n ~ 1 
and averaging over the forcing statistics. The active scalar moments obey the equation 
dt(a 2n ) = n(2n — l)F Q (a 2n ~ 2 ) (odd moments vanish by symmetry) whose solutions are 
the Gaussian moments: (a 2n ) = {2n — l)!!(F /f:) n . An equivalent derivation can be 
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Figure 9. Pdf's of active (empty boxes) and passive (full boxes) scalar fields 
normalized by their standard deviation. The active scalar pdf is indistinguishable 
from a Gaussian (solid line). 



obtained in Lagrangian terms. Following the same steps which lead from Eq. ()30j) to 
fl3*T|) it is easy to derive the following expression 

a 2n {x, t) = 2nJ&s x J&y 1 P{y l , Sl \x, t)f a { yi , s^a^fa, s x ) , (39) 
which after integrating over x and averaging over the forcing statistics reduces to 

(a 2n )(t) = n(2n-l)F f ds(a 2n " 2 )(s) , (40) 

Jo 

which, unraveling the hierarchy, yields Gaussian moments written above. In passing 
from Eq. (|39p to ([40)) we used the property that J dxP(y, s\x, t) = 1 (which is ensured 
by Eq. (fTT)|) ). and Gaussian integration by parts. 

4-4- 2. Multi-point statistics and the absence of anomalous scaling in the inverse cascade 
The results about two-point statistics of the magnetic potential and passive scalar fields 
is summarized in Fig. E3 In the inverse cascade scaling range r > if the rescaled pdf 
of S r a at different values of the separation r collapse onto the same curve indicating 
absence of anomalous scaling (Fig. ITUkh By contrast, in the scaling range r < if the 
rescaled pdfs of passive scalar increments at different r do not collapse (Fig. 110b). i.e. 
we have anomalous scaling, as expected. 

Let us now consider in quantitative terms the scaling behaviors. Since in our 
simulations the scaling range at r < if (k > kf) is poorly resolved, we do not 
enter the much debated issue of the scaling of the magnetic and velocity fields in 
this range of scales (see, e.g., [H3 EH1 H3 EBl EH1 ED])- We just mention that current 
opinions are divided between the Kolmogorov scaling 5 r v ~ r 1//3 , and the Iroshnikov- 
Kraichnan 5 r v ~ r 1//4 , corresponding to spectral behaviors such as E v {k) ~ A;~ 5 / 3 and 
E v {k) ~ A;" 3 / 2 , respectively. Both theories agree on the smooth scaling behavior for 
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5 r a/<(S r a)V /2 , 5 r c/<(5 r c)V' 2 5 r a/<(5 r a) V 2 , 5 r c/<(8 r c) 2 > 1 ' 2 

Figure 10. (a) Rescaled pdf of the normalized active scalar increments, 
<5 r a/((5 r a) 2 ) 1 / 2 , for three separations in the r > if range: r = 2£f (x) ,3£f (o), 
Atf (□). The normalized passive scalar increments pdf at the same scales (dashed 
line) and the Gaussian (solid line) are shown for comparison . In the inset the rescaled 
pdf of 6 r a for r = 3£f (solid line) is shown in the interval [—2.5 : 2.5] in linear scale, 
to emphasize the deviation from a Gaussian (dashed line). The flatness is f=a 2.76, 
significantly smaller than the Gaussian value 3. (b) Rescaled pdf of the passive scalar 
increments, S r c/ ((Src) 2 } 1 ^ 2 , for two separations in the r < It range: r = 0.3£/ (full 
triangles) and 0.6£/ (□). The rescaled pdfs of active scalar increments at the same 
scales (solid line) and a Gaussian (dashed line) are reported for comparison. 



the magnetic potential, 5 r a ~ r, that is observed numerically (see Fig. |3J). In the 
range of scales r > if (k < kf) standard dimensional arguments predict E a {k) ~ k~ 7 ^ 3 
113 EH] which is different from our findings E a {k) ~ k~ 2 (Fig. [TTJ). In real space 
this means that 5 r a ~ r 1 / 2 , which is dimensionally compatible with scaling behavior 
5 r v = [v(x + r,t) — v(x, t)] • r ~ r° for r > if (as suggested by the velocity spectrum, 
see Fig. 0}, and the Yaglom relation (5 r v(S r a) 2 ) ~ F^r jl]. As a side remark, note 
that the argument for E a (k) ~ k~ 7 ^ 3 rests on the assumption of locality for velocity 
increments, a hypothesis incompatible with the observed velocity spectrum at k < kf 
(see Fig. IU). 

It is worth remarking that the increments of active scalar d~ r a = a(x + r,t) — a(x, t) 
eventually reach a stationary state in spite of the growth of a 2 . The distribution of 
of S r a is subgaussian (Fig. [TUJ). This behaviour can be explained recalling that 5 T a is 
the difference of two Gaussian variables, which are however strongly correlated (indeed 
the main contributions will come from x + r and x inside the same island, see Fig. EJ). 
This correlation leads to cancellations resulting in the observed subgaussian probability 
density function. 

Let us now focus on the most interesting aspect, namely the absence of anomalous 
scaling in the inverse cascade range. The absence of intermittency seems to be a common 
feature of inverse cascading systems as passive scalars in compressible flows |SH ESI 
and the velocity field in two dimensional Navier Stokes turbulence [42]. This leads to 
the conjecture that an universal mechanism may be responsible for the self similarity 
of inverse cascading systems. While 2d Navier-Stokes turbulence is still far to be 
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Figure 11. Compensated power spectra of active scalar variances E a (k)k 2 at three 
different times of the evolution. Note that an increasingly clear plateau at A: < kf 
appears. 

understood, the absence of anomalous scaling in passive scalars evolving in compressible 
flows has been recently understood in terms of the collapse of the Lagrangian trajectories 
onto a unique path 23131]. Briefly, the collapse of trajectories allows to express the 
2iV-order structure function in terms of two-particle propagators, (P2)v, instead of the 
27V-particle propagator, (P2n)v While the latter may be dominated by a zero mode with 
a non trivial anomalous scaling, the former are not anomalous and lead to dimensional 
scaling. 

The above argument cannot be simply exported to the magnetic potential inverse 
cascade: first, the Lagrangian paths do not collapse onto a unique one; second, the 
correlation between the f a and v does not allow to split the averages. However, the 
former difficulty can be overcome. Indeed the property that all paths landing in the 
same point contribute the same value has important consequences also in the multipoints 
statistics. Proceeding as for the derivation of Eq. (jSH), it is possible to show that m 
trajectories are enough to calculate the product of arbitrary powers of a at m different 
points (a ni (cci,i) . . . a nm (x m , t)). In particular, structure functions S^(r) = ((S r a) N ) for 
any order N involve only two trajectories. This should be contrasted with the passive 
scalar, where the number of trajectories increases with iV and this is at the core of 
anomalous scaling of passive fields [27j . 

5. Two-dimensional Ekman turbulence 

Let us now consider the case of scalars that act on the velocity field through a functional 
dependence. Among them, probably the best known case is vorticity in two dimensions 
a = V x v. It obeys the Ekman-Navier- Stokes equation 

dta + v ■ Va = nAa — a a + f a , (41) 
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Figure 12. (a)The vorticity spectrum E a {k) ~ fc -1- '' steepens by increasing the 
Ekman coefficient a. Here a = 0.15 (+), a = 0.23 (x), a = 0.30 (©). In the inset, the 
exponent £ as a function of a. (b) Power spectra of passive scalar (x) and vorticity 
(+). Here a = 0.15. In the inset we show the same spectra compensated by /c 1+ ^ 2 . 
For details on the DNS see [T2] . 



The term —a a models the Ekman drag experienced by a thin fluid layer with the walls 
or the surrounding air |37] • 

In the absence of friction (a = 0) dimensional arguments confirmed by 

experiments |4*T] 153] and numerical simulations S42], give the following scenario. At 
large scales an inverse (non-intermittent) cascade of kinetic energy takes place with 
E v (k) ~ k~ 5 ^ 3 and E a (k) ~ k 1 ^ 3 . At small scales the enstrophy, (a 2 ), performs a 
forward cascade with E a (k) ~ k~ l and E v (k) ~ k~ 3 , meaning that the velocity field is 
smooth in this range of scales. 

In the presence of friction (a > 0) kinetic energy is removed at large scales holding 
the system in a statistically steady state and small-scale statistics is modified by the 
competition of the inertial (v ■ Va) and friction (—a a). Since these two terms have 
the same dimension due to smoothness of the velocity field, this results in nontrivial 
scaling laws for a. This effect is evident in DNS with a > ^3 E] , where the vorticity 
spectrum displays a power law steeper than in the frictionless case: E a (k) ~ 
with £ dependent on a (see Fig. IT^k). For < £ < 2 the exponent £ coincides with 
the scaling exponent, Q, of the second order structure function S%(r). Additionally, 
high-order structure functions at fixed a show anomalous scaling, ^ N£/2 [PTj . 
The spectral steepening and the presence of intermittency are observed [TT)| IT2*] also in 
passive scalars evolving in smooth flows according to 

d t c + v ■ Vc = kAc - a c + f c (42) 

Physically Eq. (pf2~j) describes the evolution of a decaying passive substance (e.g. a 
radioactive marker) 134*] . For a = the dimensional expectation E c (k) ~ k^ 1 has 
been verified in experiments [HE]- For positive a passive scalar spectra become steeper 
than k~ l and, at high wavenumbers, have the same slope of the vorticity spectrum 
(Fig. 112b). Additionally, the pdfs of passive and vorticity increments for separations 
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Figure 13. Probability density functions of vorticity differences (+) and of passive 
scalar ones (x), normalized by their respective standard deviation, at different scales 
r within the scaling range. 



inside the inertial range collapse one onto the others (Fig.ll3p. signalling that the scaling 
exponents coincide: (n = (n- m summary, the vorticity field and the passive scalar share 
the same statistical scaling properties ^TJEEE], similarly to the 2D convection case. 
However, differences may appear for odd order moments [q"o] . 

It is now interesting to understand how the equivalence of active and passive 
statistics is realized: we shall see that in this case the smoothness of the flow is a 
crucial ingredient. 

Let us start with the decaying passive scalar (}4*2*j) . First of all, it should be 
noted that the presence of non-zero friction (a > 0) regularizes the field, and there 
is absence of dissipative anomaly [Slj , even if the mechanism is different from the MHD 
one. As a consequence, in Eq. ()42)1 we can put k = and solve it by the method of 
characteristics (see Sect. I2.2|) . i.e. c(x,t) = f_ dsf c (p(s; x, t), s) exp(— a(t — s)), where 
now the path p(s;x,t),s) is unique due to the smoothness of the velocity field (note 
that E v {k) ~ k~ 3 ~^ is always steeper than k~ 3 , see also Fig. fT^k) . The integral extends 
up to — oo where the initial conditions are set. 

Passive scalar increments 5 r c = c(xi,t) — c(x2,t) (r = \x\ — a?2 1 ) , which are the 
objects we are interested in, are associated to particles pairs 

5 r c =j dse- ait - s) {f c (p( S] x u t), s)-f c (p(s; x 2 ,t), s)] . (43) 

J — oo 

It should be now noted that the integrand f c (p(s; x±,t), s)—f c (p(s; x 2 ,t), s) stays small 
as long as the separation between the two paths remains below the forcing correlation 
scale, if] while for separations larger than if it is approximately equal to a Gaussian 
random variable, X. In the latter statement we used the independence between f c and 
the particle trajectories, ensured by the passive nature of c. Therefore, Eq. (}4*3*j) can be 
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Figure 14. The scaling exponents of the passive scalar £° (+), which have been 
computed also for noninteger moments, indeed Eq. H46JI holds in general. We also 
show the exponents obtained from the separation times statistics (0) according to 
(exp[— apTe f (r)]} ~ r^p with average over about 2 x 10 5 couple of Lagrangian particles. 
The errorbars are estimated by the r.m.s. fluctuation of the local slope. In the inset 
we plot the Cramer function 5(7) computed from finite time Lyapunov exponents 
(symbols) and exit time statistics (line). 



approximated as 

5 r c^X I J ' ds e~ a{t - s) ~ Xe" x7i f {r) , (44) 



t-T lf (r) 



where %> f (r) is the time necessary for the particles pair to go from r to if backward in 
time. It is now clear that large fluctuations are associated to fast separating couples, 
1it(f) <C (% f (r)), and small fluctuations to slow ones. Moreover, since v is smooth, 
two-dimensional and incompressible pairs separation is exponential and its statistics 
is described by a single finite-time Lyapunov exponent j2Z], 7- It is related to the 
separation time through the relation: 

£ f = re lTe f {r) . (45) 

Large deviations theory states that at large times the random variable 7 is distributed as 
P(7, i) ~ t 1 ! 2 exp[— S{^i)t\. 5(7), the so-called Cramer function [HE], is positive, concave 
and has a quadratic minimum, 5(A) = in A, the maximum Lyapunov exponent. From 
(I44j) and (|43|) along with with the expression for P( / y,t), structure functions can be 
computed as 

AMD " »(£)", (46) 

where = min 7 {A r , [Na + S^)]/^}. Fig. HU shows that this prediction is verified 
by numerical simulations. The anomalous exponents can therefore be evaluated in 
terms of the flow properties through the Cramer funtion S^). 
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Since = one ma y be tempted to apply straightforwardly to vorticity the 
same argument used for the passive scalar. However, the crucial assumption used in 
the derivation of Eq. (jUjj) is the statistical independence of particle trajectories and the 
forcing. That allows to consider % f (r) and X as independent random variables. While 
this is obviously true for c, it is clear that the vorticity forcing, f a , may influence the 
velocity statistics and, as a consequence, the particle paths. In other terms, for the 
vorticity field we cannot a priori consider X and % f (r) as independent. 

Nevertheless, it is possible to argue in favour of the validity of Eq. (jlHjl also for 
vorticity. Indeed, the random variable % f (r) is essentially determined by the evolution 
of the strain along the Lagrangian paths for times t — % f (r) < s < t, whilst X results 
from the scalar input accumulated at times s < t — % f (r). Since the strain correlation 
time is basically a -1 , for TtAr) ^> a~ l it is reasonable to assume that % f {r) and X are 
statistically independent. Eq. (pIK]) allows to translate the condition % } {r) a^ 1 into 
r <C £jexp(— j/a), and to conlcude that at small scales we can safely assume that the 
vorticity behaves as a passive field. 

We conclude this section with two remarks. First, a difference between the present 
scenario and the one arisen in two-dimensional convection should be pointed out. Here 
the scaling exponents depend on the statistics of the finite time Lyapunov exponents, 
which in turn depends on the way the vorticity is forced. As a consequence universality 
may be lost. Second, the smoothness of the velocity field plays a central role in the 
argument used to justify the equivalence of the statistics of a and c in this system. For 
rough velocity fields % f {r) is basically independent of r for r << if and the argument 
given in this section would not apply. 

6. Turbulence on fluid surfaces 

6.1. Surface Quasi Geostrophic turbulence 

The study of geophysical fluids is of obvious importance for the understanding of weather 
dynamics and global circulation. Taking advantage of stratification and rotation, 
controlled approximations on the Navier-Stokes equations are usually done to obtain 
more tractable models. For example, the motion of large portions of the atmosphere 
and ocean have a stable density stratification, i.e. lighter fluid sits above heavier one. 
This stable stratification, combined with the planetary rotation and the consequent 
Coriolis force, causes the most energetic motion to occur approximately in horizontal 
planes and leads to a balance of pressure gradients and inertial forces. This situation 
is mathematically described by the Quasi- Geostrophic equations pHIH]. On the lower 
boundary, the surface of the earth or the bottom of the ocean, the vertical velocity has 
to be zero and this further simplifies the equations. Assuming that the fluid is infinitely 
high and that the free surface is flat, the dynamics is then described by the Surface 
Quasi-Geostrophic equation (SQG) [HJ EH]- It governs the temporal variations of the 
density field, our active scalar a throughout this section. For an ideal fluid the density 
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variation will be proportional to the temperature variation and one can use the potential 
temperature as fundamental field. 

The density fluctuations evolve according to the transport equation 

d t a + v ■ Va = «Aa + f a , (47) 

and the velocity field is functionally related to a. In terms of the stream function ip, the 
density a is obtained as a(x, t) = ( y —A) 1 ^ 2 -ip(x, t), and v = (d y , —d x )ip. In Fourier space 
the link is particularly simple 



«(fc,t) = -i a(M ' 



(48) 



In the early days of computer simulations eqs. (|4~7jl and (}4*%j) had been used for weather 
forecasting. Recently they received renewed interest in view of their formal analogy 
with the 3d Navier-Stokes equations |H0Jini], and as a model of active scalar [9*1 ITU] 
It is instructive to generalize ("""""|) as 

v(k,t) = -i(^,-^)a(k,t). (49) 



k z k z / 

For z = 1 one recovers (J48|) . while z = 2 corresponds to the 2d Navier-Stokes equation 
and a is the vorticity (see the previous section). The value of z tunes the degree of 
locality. The case of interest here is z = 1, corresponding to the same degree of locality 
as in 3d turbulence [§]. 

Eqs. (j47|) . (j49j) have two inviscid unforced quadratic invariants E = — J dxaip and 
Q = J dxa 2 , which for z = 2 correspond to the energy and the enstrophy, respectively. 
Notice that, for z = 1, the fields v and a have the same dimension and E a (k) = E v (k). 
Quasi-equilibrium arguments predict an inverse cascade of E and a direct cascade of 
Q. If kf identifies the injection wave number dimensional arguments give the expected 
spectral behavior [""]: 



7 , 4z 
" 3" 1 " 3 



For z = 2 one recovers the 2d Navier-Stokes expectation [SH Here, we are 

interested in the range k > kf and z — 1, so that E a {k) = E v (k) ~ k~ 5 ^ 3 as 
in 3d turbulence. Assuming absence of intermittency, the spectrum E v {k) ~ k~ 5 ^ 3 
would imply 5 r v ~ r 1 / 3 . Therefore, for a passive scalar one expects the standard 
phenomenology: an intermittent direct cascade with E c {k) ~ k~ 5 ^ 3 to the Oboukov- 
Corrsin-Kolmogorov type of arguments [3]. 

6.1.1. Direct numerical simulations settings We performed a set of high resolution 
direct numerical simulations of Eqs. ()47|) and (}4*%|) along with the passive scalar equation 
(J2J. Integration has been performed by means of a 2/3 de-aliased pseudo-spectral 
method on a 2n x 27r doubly periodic square domain with Af 2 collocation points. 
Time integration has been done with a second-order Adam-Bashforth or Runge-Kutta 
algorithms, approriately modified to exactly integrate the dissipative terms. The latter 
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Table 1. Summary of settings in DNS. 



run 


N 2 


P 


'7 


F 


*/ 


ja 


4096 2 


2 




B 


44-48 


2 b 


2048 2 


2 


1/2 


B 


2-6 


;i<> 


2048 2 


2 





G 


5 




2048 2 


2 


1/2 


G 


■5 


5 C 


2048 2 


2 





NG 


5 



a Runs 1,2 are forced according to Eq. (|23[1 . and run 1 is without any friction term. 
b Run 3,4 are forced according to Eq. I|27|) . 

c Run 5 is forced with the non-Gaussian forcing discussed in Sect. 16.31 Lower (512 2 ) 
resolutions runs with several settings both for the dissipative and friction terms have 
been also performed 

ones, as customary, have been replaced by a hyper-diffusive term that in Fourier space 
reads —k 2p a(k,t). Since the system is not stationary due to the inverse cascade of E, 
we added an energy sink at large scales in the form —k~ q d(k, t). In order to evaluate its 
possible effect on inertial quantities, a very high resolution (4096 2 ) DNS has been done 
without any energy sink at large scales. 

A summary of the numerical settings can be found in table We considered 
different scalar inputs: (G) a Gaussian ^-correlated in time forcing as (J2HJ); (B) a 
5-correlated in time one restricted to a few wavenumbers shells as (J2ZJ); (NG) a 
non-Gaussian non 5-correlated in time one suited to produce non zero three points 
correlations for the scalar fields (see below in Sect. EEHJ). 

6.2. Active and passive scalar statistics in SQG turbulence 

Let us start looking at the measured scalar spectra in order to check the dimensional 
predictions. In Fig. we summarize the spectral behavior for a and c in two different 
simulations (i.e. runs 1 and 3). 

Two observations are in order: 

• The scalar spectra E a {k) is steeper than the dimensionally expected k~ 5 ^ 3 . The 
slope does not seem to depend on the injection mechanism. 

• The active and passive scalars, both performing a direct cascade, are different 
already at the level of the spectral slopes. In particular, c is much rougher than a. 

The differences between a and c are evident also from Fig. where two 
simultaneous snapshots of the fields are displayed. A direct inspection of the fields 
confirms that c is more rough than a, and resembles a passive scalar in a smooth flow. 
Moreover, in a we observe the presence of large coherent structures at the scale of the 
forcing. These are actually long-lived, slowly-evolving structures that strongly impair 
the convergence of the statistics for high-order structure functions. In the following 
we thus limit ourselves to a comparison of the low-order statistics of a and c. That 
is however sufficient to appreciate the differences beween active and passive scalar 
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Figure 15. (a) Power spectra of active (red) and passive (blue) scalar variances 
E a (k) — irk\a(k,t)\ 2 and E c (k) — irk\c(k, t)\ 2 for run 1, the inset shows the active, 
and passive energy fluxes IT a)C (fc). The dashed lines corresponds to the best fitted 
spectral slopes E a (k) ~ ^ L8±al and E c (k) ~ fc-i.i5±o.05_ Note that thc fitted 
slope for E a (k) is better recovered at higher wavenumbers, close to the energy peak 
a steeper (~ k~ ) is observed, (b) The same but for run 3. The left inset shows 
the scalar fluxes. The dashed lines indicates the spectral slopes E a (k) ~ fc-i-8±0.i anc j 
E c {k) ~ fc-i-i7±0.05_ jj ere th e scaling range is wider than in run 1. For both runs 
active and passive spectra have been shifted for visualization purposes. In thc other 
runs we observe the same qualitative and quantitative features, within the error bars. 
In particular, run 1 and run 2 give almost indistinguishable spectral slopes, meaning 
that the large scale energy sink does not enter too much into the inertial range, as 
expected since the velocity field is rough at small scales. Run 5 produced spectral 
slopes indistinguishable from those of runs 3 and 4. However, it should be noted that 
run 1 and run 2 display a steeper spectra close to the forcing wave numbers, this 
feature is absent in the other runs and may be due to finite size effects, sec text and 
also Fig. El 



statistics. 

The deviation from the dimensional expectation for the spectral slope of the active 
field was already observed in previous numerical studies and it is possibly due to 
intermittency in the active scalar and velocity statistics. Indeed, the rescaled pdfs of 
the increments do not collapse (not shown). Concerning the universality of anomalous 
exponents with respect to the forcing statistics, we observe that the spectral slopes do 
not depend sensitively on the forcing statistics. However, the forcing (B) (see table GJ 
induces a steeper spectrum close to the forcing scale, while a universal slope is apparently 
recovered at large wavenumbers. (Fig. IT5k). The scaling of S% (r) seems to be more 
sensitive to the forcing (see Fig. IT7]) . These discrepancies may be due to finite size 
effects which are more severe for forcing (B) which is characterized by a slower decay 
of the spatial correlations. Unfortunately, it is difficult to extract reliable information 
on the high-order statistics. Indeed the presence of slowly evolving coherent structures 
induces a poorly converging statistics for high-order structure functions. Therefore, we 
cannot rule out the possibility of a dependence of active scalar exponents on forcing 
statistics. 
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Figure 16. Snapshots of the active (left) and passive (right) scalar (at 512 2 resolution). 
Note the presence of coherent structures in the active field, which are slow evolving. 
The passive field displays filamental like features indicating that it is more rough than 
the active one. 



A result that emerges beyond any doubt is that active and passive scalars behave 
differently, as shown in Figs. E3 EI] and El This is confirmed by the differences in the 
pdfs of 5 r a and 5 r c at various scales within the inertial range (Fig.UHjl. The single-point 
pdfs of a and c are different as well (not shown). 

It is worth noticing that the behavior of the passive scalar deviates from naive 
expectations. We observe E c (k) ~ k~ 1A5 whereas a dimensional argument based on the 
observed velocity spectrum (E v (k) = E a {k) ~ k~ 18 ) and on the Yaglom relation [I] 




Figure 17. (a) Second order structure function for the active (empty boxes) and the 
passive (full circles) scalar in run 2. The straight line shows the scaling ~ r 1 for the 
S%(r) (the same exponent is observed in run 1). The passive scalar does not display 
a neat scaling behavior. In the inset we show the correlation function for the passive 
scalar, C|(r), in log-lin scale, (see text), (b) The same but for run 3. Here the straight 
dashed line indicates the slope r 08 . S%(r) in run 5 has a slope ~ r °- 84 ; compatible 
with the ones observed in runs 3 and 4 within statistical errors. 
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Figure 18. Pdf of active (red) and passive (blue) scalar differences normalized by 
their respective standard deviations for four different scales within the inertial range. 
The dotted lines are Gaussian pdfs drawn for comparison. Data refer to run 3, the 
other runs show similar results. 



((S r v (S r c) 2 ) ~ r) would give E a (k) ~ A;^ 16 . The violation of the dimensional prediction 
is even more striking looking at S%(r) in Fig. El This feature is reminiscent of some 
experimental investigations of passive scalars in turbulent flows, see e.g. Refs. jH2 EH] 
and references therein, where shallow scalar spectra are observed for the scalar even 
at those scales where the velocity field displays a K41 spectrum. It is likely that the 
presence of coherent structures (see Fig.lT^j) leads to persistent regions where the velocity 
field has a smooth (shear-like) behavior. This suggest a two-fluid picture: a slowly 
evolving shear-like flow, with superimposed faster turbulent fluctuations. Under those 
conditions, one may expect that the particle pairs separate faster than expected, leading 
to shallower passive scalar spectra [S3]. 

In conclusion, even though both passive and active scalar perform a direct cascade, 
their statistical properties are definitely different. 

6.3. Scaling and geometry 

Let us now compare the two scalars field by investigating the three-point correlation 
functions of the active field, Q (a?i, x 2 , x 3 ) = (a(xi,t)a(x 2 , t)a(x 3 , t)), and of the passive 
one, C 3 (xi,x 2 ,x 3 ) = (c(xx, t)c(x 3 , t)c(x 3 , t)). This allows to compare the scaling 
properties of the correlation function by measuring the dependence on the global size 
of the triangle identified by the three points, R 2 = Ylii<j x % (being Xij = \x{ — Xj\). We 
will also investigate the geometrical dependence of C 3 . 

It has to be noted that for 5-correlated Gaussian forcings as (J23)l and (J27j) . C 3 and 
C 3 are identically zero. Therefore we need a different forcing statistics in order to study 
three-point correlations. A possibility is to break the rotational symmetry of the system 
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by an anisotropic forcing, e.g., a mean gradient (J24)) as in |28j . However, that choice 
leads inevitably to a = c: the equations are identical for a and c so that the difference 
field a — c will decay out. We then choose to force the system as follows. The two inputs 
f a and f c are 

f atC (x, t) = gl c (x, t) - J dygl ;C (y, t) (51) 
where g a c is a homogeneous, isotropic and Gaussian random field with correlation 



(g i (x,t)g j (x',t')) = 5 i!j g(\x-x l \/£ f )e ^ , (52) 

where i,j = a,c, if is the forcing scale, Ty the forcing correlation time and G{r) oc 
Go exp(— r 2 /2). The time correlation is imposed by performing an independent Ornstein- 
Uhlenbeck process at each Fourier mode, i.e. integrating the stochastic differential 
equation dcji(k,t) = —l/Tfcji(k,t)dt + y^2Godt/rfdwi(k, t) (where duii are zero mean 
Gaussian variables with (dwi(k,t)dvjj(k' ,t')) = 5 it j5(t — t')5{k — k')). If Tf <C %. by the 
central limit theorem a Gaussian statistics is recovered. Therefore, we fixed Tf ~ 0(% f ) 
in our DNS. The advantages of this choice are that it preserves the isotropy and gives 
analytical control on the forcing correlation functions. 

Let us now see how the triangle identified by the three points, x = (xi,x 2 ,x 3 ) 
can be parameterized. In two dimensions we need 3d = 6 variables to define a 
triangle. Since the correlation functions should possess all the statistical symmetries 
of the system the number of degrees of freedom is reduced. In particular, translational 
invariance ensures no dependence on the position of the center of mass of the triangle, 
(xi + x 2 + x 3 )/?>. The correlation function is thus a function of the separation vectors 
among the 3 points, i.e. C^' c (x) = C^ ,c (xi 2 , x 23 , x 3 i). Additionally, isotropy implies that 
a rigid rotation of the triangle has no effect on the value of C 3 , so that three variables 
suffice: the global size of the triangle R, and two parameters that define its shape. In 
terms of the Euler parametrization jHSHEElEZj, upon defining p\ = (x\ — a? 2 )/\/2 and 
P2 — ( x i + x 2 — 2a; 3 ) / the shape of the triangle is given in terms of the two variables: 

2 pi ■ p 2 



2pi x p 2 1 - 

w = y — - tan 

R 2 A 2 



Px-pl 



(53) 



where \w\ is the ratio of the area of the triangle divided by the area of the equilateral 
triangle having the same size, R 2 = p\ + p\. x has not a simple geometrical 
interpretation. Some shapes corresponding to a few (w, x) ar e shown in Fig. fTTH 
The three-point correlation function can be decomposed as [SB I2Z| 

where TZ,T are the reducible and irreducible components, respectively. This means 
that C3 can be expressed as the sum of a part that depends on three points (hereafter 
denoted as C3) and a part that depends on two points (hereafter denoted C 2 ), i.e. 
C 3 (:ri, x 2 , x 3 ) = C^(xi, x 2 , x 3 ) + C 2 (x 1 , x 2 ) + C 2 (x 2l x 3 ) + C 2 (x 3 , xi). The reducible, 
1Z{Cs} = C 2 (xi, x 2 )+C 2 (x 2 , x 3 )+C 2 (x 3 , Xi) and irriducible part T{C 3 } = 6*3(^1, #2, #3) 
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Figure 19. Contour plots in the \ ~ w plane of the three-point shape function, 
for passive (left) and active (right) scalar. 



are in general characterized by different scaling properties as a function of the triangle 
size R, and different geometrical (shape) dependencies (53 EH] • 

In terms of the variables R, w and x, the reducible and irreducible components take 
the following form 

T{Cr(x)} =hl c (R)<S>l c ( X ,w), (55) 

where the function h(R) is expected to have a scaling dependence in the inertial range. 

For the passive scalar, the scaling behavior of the reducible part, h^-(R), is 
dominated by the dimensional scaling imposed by the balance with the forcing, while 
the scaling of the irreducible part, ti^(R), is given by the zero modes [31J. With 
a finite-correlated pumping the forcing statistics may in principle contribute to the 
irreducible part j22j; however in our case these hypothetical contributions seem small, 
if not absent. For the active scalar it is not possible to make any a priori argument 
to predict the scaling behavior of the reducible and irreducible part. As far as the 
geometrical dependence is concerned, <&^> J and , it is very difficult to separate 
the two contributions. In our case, in agreement with the results obtained for the 
Kraichanan model [BE], the reducible part turns out to be the leading contribution, 
so that we only present the shape dependence of the full correlation functions which 
basically coincides with the reducible one. Hereafter we will then use $ ajC dropping the 
indices which distinguish the two contributions. 

Let us start by studying the shape dependence of the correlation function for the 
active, $ a (x w )i an d passive scalars, $ c (x, w). Exploiting the invariance under arbitrary 
permutations of the three veritices of the triangle we can reduce the configuration 
space to x £ [ — 7r/6 : 7r/6], and w £ [0 : 1] going from degenerate (collinear points) 
to equilateral triangles. The function $>(x,w) is antiperiodic in x with period 7r/3 
H23 EZ| • In Fig- EH the functions $ c and $ c are shown. They have been measured 
for a fixed size R within the inertial range. The two functions display similar qualitative 
features: the intensity grows going from equilateral to collinear triangles, and the 
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Figure 20. (a) h n {R) vs R for the active (o) and passive (□) scalar. The dashed 
lines are the second order structure functions, S%' c , multiplied by a factor to allow for 
a comparison. The solid line gives the fitted slope, jj 1 - 23±0 - 05 ) for ti^{R). The slope 
of S% is ~ r -84±o.05 an( j j g conl p a tible with those obtained with the forcing (G), see 
Fig. I17fl . The passive curves have been shifted for visualization purposes, (b) C (R) 
vs R for the active (o) and passive (□) scalar. For c the scaling R 1 - 7 ^ - 2 (solid line) is 
measured. The signal for the active scalar is very low and no scaling behavior can be 
detected. The statistics has been computed averaging over about 65 frames separated 
by half eddy turnover time. 



maximum is realized for almost degenerate triangles (%, w) = (0,0). $ c is invariant 
for x ~> ~Xj which corresponds to reflection with respect to an axis. This symmetry is 
a consequence of the equation of motion and the fact that f c — > f c under this symmetry 
transformation. The active scalar displays a weak breaking of this symmetry, as a 
consequence of the fact that it is a pseudoscalar while f a is not. 

Let us now turn our attention to the scaling behavior. Since the reducible scaling 
behavior is the leading one it is simply obtained by fixing a certain shape for the triangle 
(we did for several choices of (%, w)) and varying its size. The results are shown in in 
Fig. nOh . The observation is that h£(R) and ti^(R) are different. Moreover, while 
h™(R) is practically parallel to S%(R), for the active scalar we observe that h™(R) is 
not scaling as S%(R). 

The measure of the subleading, irreducible part is more involved, and we proceed 
as follows. We fix a reference triangle shape (x, w) and set the origin in the center 
of mass of the triangle. Now a;, indicates the position of the vertex % of the 
triangle. We define d\ = di(\) (A > 1) as the dilation operator which transforms 
the triangle (xi,X2,Xs) in (Acci, x 2 , X3). Analogous definitions hold for the other 
vertices. Obviously <ij(l) = / is the identity. We then consider the composite operator 
T){\) = ^1(^2^3 — did 2 — d 2 d3 — d^di + d\ + d 2 + (I3 — /. By direct substitution it is easily 
seen that in Z>(A)C 3 (£Ci, x 2 , x 3 ) all the reducible terms are cancelled, and only a linear 
superposition of terms involving the irreducible parts survives Therefore, the average of 
V(\)c(xi)c(x2)c(x 3 ) for triangles of different sizes R but of fixed shape w) will give 
the scaling of the irreducible part of the three-point function. ([55)1 . A similar procedure 
has been done for a. We used as reference configuration a collinear triangle with two 
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degenerate vertices, i.e. (x, w) = (0, 0). This choice reduces the number of cancellations. 
Additionally, this configuration corresponds to the region of stronger gradients in the 
function $ (see Fig. HH|) . yielding a higher signal to noise ratio. By varying A we tested 
the robustness of the measured scaling. 

In Fig. EDI we present the results on the scale dependence of h^(R) and h^(R). 
Clearly h^(R) has too low a signal to identify any scaling behavior, while the passive 
scalar scales fairly well and we measured h^(R) ~ R 17 , confirming that the irreducible 
part is subleading. 

A couple of final remarks are in order. First, the fact that K^(R) does not scale 
as S*2 (R) is the signature of the correlations between f a and the particle propagator. 
Indeed, for the passive scalar ti^(R) oc S^R) is a straightforward consequence of the 
independence of v and f c . Second, h^(R) scales differently from h^(R): this rules out 
the possibility of establishing simple relationships between active and passive scalar 
statistics (see [T3] for a related discussion). 

7. Conclusions and perspectives 

In summary, we have investigated the statistics of active and passive scalars transported 
by the same turbulent flow. We put the focus on the issue of universality and scaling. 
In this respect the passive scalar problem is essentially understood. Conversely, the 
active case is by far and away a challenging open problem. The basic property 
that make passive scalar turbulence substantially simpler is the absence of statistical 
correlations between scalar forcing and carrier flow. On the contrary, the hallmark 
of active scalars is the functional dependence of velocity on the scalar field and thus 
on active scalar pumping. Yet, when correlations are sufficiently weak, the active 
scalar behaves similarly to the passive one: this is the case of two-dimensional thermal 
convection and Ekman turbulence. However, this appears to be a nongeneric situation, 
and the equivalence between passive and active scalar is rooted in special properties of 
those systems. Indeed, different systems as two-dimensional magnetohydrodynamics or 
surface quasi-geostrophic turbulence are characterized by a marked difference between 
passive and active scalar statistics. This poses the problem of universality in active 
scalar turbulence: if the forcing is capable of influencing the velocity dynamics, how can 
scaling exponents be universal with respect to the details of the injection mechanism 
? So far, a satisfactory answer is missing. As of today, numerical experiments favor 
the hypothesis that universality is not lost in a number of active scalars, but further 
investigation is needed to elucidate this fundamental issue. 
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